3.1.50 \(\int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx\) [50]

3.1.50.1 Optimal result
3.1.50.2 Mathematica [A] (verified)
3.1.50.3 Rubi [A] (verified)
3.1.50.4 Maple [B] (warning: unable to verify)
3.1.50.5 Fricas [F(-1)]
3.1.50.6 Sympy [F]
3.1.50.7 Maxima [F]
3.1.50.8 Giac [F]
3.1.50.9 Mupad [F(-1)]

3.1.50.1 Optimal result

Integrand size = 35, antiderivative size = 398 \[ \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx=-\frac {2 \sqrt {c+d} \operatorname {EllipticPi}\left (\frac {a (c+d)}{(a+b) c},\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right ),\frac {(a-b) (c+d)}{(a+b) (c-d)}\right ) \sec (e+f x) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (1+\sin (e+f x))}{(c-d) (a+b \sin (e+f x))}} (a+b \sin (e+f x))}{a \sqrt {a+b} c f}-\frac {2 b \sqrt {a+b} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right ) \sec (e+f x) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (1+\sin (e+f x))}{(a-b) (c+d \sin (e+f x))}} (c+d \sin (e+f x))}{a \sqrt {c+d} (b c-a d) f} \]

output
-2*EllipticPi((a+b)^(1/2)*(c+d*sin(f*x+e))^(1/2)/(c+d)^(1/2)/(a+b*sin(f*x+ 
e))^(1/2),a*(c+d)/(a+b)/c,((a-b)*(c+d)/(a+b)/(c-d))^(1/2))*sec(f*x+e)*(a+b 
*sin(f*x+e))*(c+d)^(1/2)*(-(-a*d+b*c)*(1-sin(f*x+e))/(c+d)/(a+b*sin(f*x+e) 
))^(1/2)*((-a*d+b*c)*(1+sin(f*x+e))/(c-d)/(a+b*sin(f*x+e)))^(1/2)/a/c/f/(a 
+b)^(1/2)-2*b*EllipticF((c+d)^(1/2)*(a+b*sin(f*x+e))^(1/2)/(a+b)^(1/2)/(c+ 
d*sin(f*x+e))^(1/2),((a+b)*(c-d)/(a-b)/(c+d))^(1/2))*sec(f*x+e)*(c+d*sin(f 
*x+e))*(a+b)^(1/2)*((-a*d+b*c)*(1-sin(f*x+e))/(a+b)/(c+d*sin(f*x+e)))^(1/2 
)*(-(-a*d+b*c)*(1+sin(f*x+e))/(a-b)/(c+d*sin(f*x+e)))^(1/2)/a/(-a*d+b*c)/f 
/(c+d)^(1/2)
 
3.1.50.2 Mathematica [A] (verified)

Time = 1.62 (sec) , antiderivative size = 374, normalized size of antiderivative = 0.94 \[ \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx=\frac {2 \sec (e+f x) \left (-\frac {(c+d) \operatorname {EllipticPi}\left (\frac {a (c+d)}{(a+b) c},\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right ),\frac {(a-b) (c+d)}{(a+b) (c-d)}\right ) \sqrt {\frac {(b c-a d) (-1+\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (1+\sin (e+f x))}{(c-d) (a+b \sin (e+f x))}} (a+b \sin (e+f x))}{c}-\frac {b (a+b) \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right ) \sqrt {\frac {(-b c+a d) (-1+\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {\frac {(-b c+a d) (1+\sin (e+f x))}{(a-b) (c+d \sin (e+f x))}} (c+d \sin (e+f x))}{b c-a d}\right )}{a \sqrt {a+b} \sqrt {c+d} f} \]

input
Integrate[Csc[e + f*x]/(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]]) 
,x]
 
output
(2*Sec[e + f*x]*(-(((c + d)*EllipticPi[(a*(c + d))/((a + b)*c), ArcSin[(Sq 
rt[a + b]*Sqrt[c + d*Sin[e + f*x]])/(Sqrt[c + d]*Sqrt[a + b*Sin[e + f*x]]) 
], ((a - b)*(c + d))/((a + b)*(c - d))]*Sqrt[((b*c - a*d)*(-1 + Sin[e + f* 
x]))/((c + d)*(a + b*Sin[e + f*x]))]*Sqrt[((b*c - a*d)*(1 + Sin[e + f*x])) 
/((c - d)*(a + b*Sin[e + f*x]))]*(a + b*Sin[e + f*x]))/c) - (b*(a + b)*Ell 
ipticF[ArcSin[(Sqrt[c + d]*Sqrt[a + b*Sin[e + f*x]])/(Sqrt[a + b]*Sqrt[c + 
 d*Sin[e + f*x]])], ((a + b)*(c - d))/((a - b)*(c + d))]*Sqrt[((-(b*c) + a 
*d)*(-1 + Sin[e + f*x]))/((a + b)*(c + d*Sin[e + f*x]))]*Sqrt[((-(b*c) + a 
*d)*(1 + Sin[e + f*x]))/((a - b)*(c + d*Sin[e + f*x]))]*(c + d*Sin[e + f*x 
]))/(b*c - a*d)))/(a*Sqrt[a + b]*Sqrt[c + d]*f)
 
3.1.50.3 Rubi [A] (verified)

Time = 1.10 (sec) , antiderivative size = 398, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.143, Rules used = {3042, 3426, 3042, 3297, 3424}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {1}{\sin (e+f x) \sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}dx\)

\(\Big \downarrow \) 3426

\(\displaystyle \frac {\int \frac {\csc (e+f x) \sqrt {a+b \sin (e+f x)}}{\sqrt {c+d \sin (e+f x)}}dx}{a}-\frac {b \int \frac {1}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}dx}{a}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\int \frac {\sqrt {a+b \sin (e+f x)}}{\sin (e+f x) \sqrt {c+d \sin (e+f x)}}dx}{a}-\frac {b \int \frac {1}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}dx}{a}\)

\(\Big \downarrow \) 3297

\(\displaystyle \frac {\int \frac {\sqrt {a+b \sin (e+f x)}}{\sin (e+f x) \sqrt {c+d \sin (e+f x)}}dx}{a}-\frac {2 b \sqrt {a+b} \sec (e+f x) (c+d \sin (e+f x)) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (\sin (e+f x)+1)}{(a-b) (c+d \sin (e+f x))}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right )}{a f \sqrt {c+d} (b c-a d)}\)

\(\Big \downarrow \) 3424

\(\displaystyle -\frac {2 b \sqrt {a+b} \sec (e+f x) (c+d \sin (e+f x)) \sqrt {\frac {(b c-a d) (1-\sin (e+f x))}{(a+b) (c+d \sin (e+f x))}} \sqrt {-\frac {(b c-a d) (\sin (e+f x)+1)}{(a-b) (c+d \sin (e+f x))}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}{\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}\right ),\frac {(a+b) (c-d)}{(a-b) (c+d)}\right )}{a f \sqrt {c+d} (b c-a d)}-\frac {2 \sqrt {c+d} \sec (e+f x) (a+b \sin (e+f x)) \sqrt {-\frac {(b c-a d) (1-\sin (e+f x))}{(c+d) (a+b \sin (e+f x))}} \sqrt {\frac {(b c-a d) (\sin (e+f x)+1)}{(c-d) (a+b \sin (e+f x))}} \operatorname {EllipticPi}\left (\frac {a (c+d)}{(a+b) c},\arcsin \left (\frac {\sqrt {a+b} \sqrt {c+d \sin (e+f x)}}{\sqrt {c+d} \sqrt {a+b \sin (e+f x)}}\right ),\frac {(a-b) (c+d)}{(a+b) (c-d)}\right )}{a c f \sqrt {a+b}}\)

input
Int[Csc[e + f*x]/(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]]),x]
 
output
(-2*Sqrt[c + d]*EllipticPi[(a*(c + d))/((a + b)*c), ArcSin[(Sqrt[a + b]*Sq 
rt[c + d*Sin[e + f*x]])/(Sqrt[c + d]*Sqrt[a + b*Sin[e + f*x]])], ((a - b)* 
(c + d))/((a + b)*(c - d))]*Sec[e + f*x]*Sqrt[-(((b*c - a*d)*(1 - Sin[e + 
f*x]))/((c + d)*(a + b*Sin[e + f*x])))]*Sqrt[((b*c - a*d)*(1 + Sin[e + f*x 
]))/((c - d)*(a + b*Sin[e + f*x]))]*(a + b*Sin[e + f*x]))/(a*Sqrt[a + b]*c 
*f) - (2*b*Sqrt[a + b]*EllipticF[ArcSin[(Sqrt[c + d]*Sqrt[a + b*Sin[e + f* 
x]])/(Sqrt[a + b]*Sqrt[c + d*Sin[e + f*x]])], ((a + b)*(c - d))/((a - b)*( 
c + d))]*Sec[e + f*x]*Sqrt[((b*c - a*d)*(1 - Sin[e + f*x]))/((a + b)*(c + 
d*Sin[e + f*x]))]*Sqrt[-(((b*c - a*d)*(1 + Sin[e + f*x]))/((a - b)*(c + d* 
Sin[e + f*x])))]*(c + d*Sin[e + f*x]))/(a*Sqrt[c + d]*(b*c - a*d)*f)
 

3.1.50.3.1 Defintions of rubi rules used

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3297
Int[1/(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(c_) + (d_.)*sin[(e_ 
.) + (f_.)*(x_)]]), x_Symbol] :> Simp[2*((c + d*Sin[e + f*x])/(f*(b*c - a*d 
)*Rt[(c + d)/(a + b), 2]*Cos[e + f*x]))*Sqrt[(b*c - a*d)*((1 - Sin[e + f*x] 
)/((a + b)*(c + d*Sin[e + f*x])))]*Sqrt[(-(b*c - a*d))*((1 + Sin[e + f*x])/ 
((a - b)*(c + d*Sin[e + f*x])))]*EllipticF[ArcSin[Rt[(c + d)/(a + b), 2]*(S 
qrt[a + b*Sin[e + f*x]]/Sqrt[c + d*Sin[e + f*x]])], (a + b)*((c - d)/((a - 
b)*(c + d)))], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && N 
eQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && PosQ[(c + d)/(a + b)]
 

rule 3424
Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]/(sin[(e_.) + (f_.)*(x_)]*Sqr 
t[(c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Simp[-2*((a + b*Sin[ 
e + f*x])/(c*f*Rt[(a + b)/(c + d), 2]*Cos[e + f*x]))*Sqrt[(-(b*c - a*d))*(( 
1 - Sin[e + f*x])/((c + d)*(a + b*Sin[e + f*x])))]*Sqrt[(b*c - a*d)*((1 + S 
in[e + f*x])/((c - d)*(a + b*Sin[e + f*x])))]*EllipticPi[a*((c + d)/(c*(a + 
 b))), ArcSin[Rt[(a + b)/(c + d), 2]*(Sqrt[c + d*Sin[e + f*x]]/Sqrt[a + b*S 
in[e + f*x]])], (a - b)*((c + d)/((a + b)*(c - d)))], x] /; FreeQ[{a, b, c, 
 d, e, f}, x] && NeQ[b*c - a*d, 0] && NeQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0 
]
 

rule 3426
Int[1/(sin[(e_.) + (f_.)*(x_)]*Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*S 
qrt[(c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> Simp[-b/a   Int[1/ 
(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]]), x], x] + Simp[1/a   In 
t[Sqrt[a + b*Sin[e + f*x]]/(Sin[e + f*x]*Sqrt[c + d*Sin[e + f*x]]), x], x] 
/; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && (NeQ[a^2 - b^2, 0] 
|| NeQ[c^2 - d^2, 0])
 
3.1.50.4 Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(24289\) vs. \(2(368)=736\).

Time = 7.12 (sec) , antiderivative size = 24290, normalized size of antiderivative = 61.03

method result size
default \(\text {Expression too large to display}\) \(24290\)

input
int(1/sin(f*x+e)/(a+b*sin(f*x+e))^(1/2)/(c+d*sin(f*x+e))^(1/2),x,method=_R 
ETURNVERBOSE)
 
output
result too large to display
 
3.1.50.5 Fricas [F(-1)]

Timed out. \[ \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx=\text {Timed out} \]

input
integrate(1/sin(f*x+e)/(a+b*sin(f*x+e))^(1/2)/(c+d*sin(f*x+e))^(1/2),x, al 
gorithm="fricas")
 
output
Timed out
 
3.1.50.6 Sympy [F]

\[ \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx=\int \frac {1}{\sqrt {a + b \sin {\left (e + f x \right )}} \sqrt {c + d \sin {\left (e + f x \right )}} \sin {\left (e + f x \right )}}\, dx \]

input
integrate(1/sin(f*x+e)/(a+b*sin(f*x+e))**(1/2)/(c+d*sin(f*x+e))**(1/2),x)
 
output
Integral(1/(sqrt(a + b*sin(e + f*x))*sqrt(c + d*sin(e + f*x))*sin(e + f*x) 
), x)
 
3.1.50.7 Maxima [F]

\[ \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx=\int { \frac {1}{\sqrt {b \sin \left (f x + e\right ) + a} \sqrt {d \sin \left (f x + e\right ) + c} \sin \left (f x + e\right )} \,d x } \]

input
integrate(1/sin(f*x+e)/(a+b*sin(f*x+e))^(1/2)/(c+d*sin(f*x+e))^(1/2),x, al 
gorithm="maxima")
 
output
integrate(1/(sqrt(b*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)*sin(f*x + e 
)), x)
 
3.1.50.8 Giac [F]

\[ \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx=\int { \frac {1}{\sqrt {b \sin \left (f x + e\right ) + a} \sqrt {d \sin \left (f x + e\right ) + c} \sin \left (f x + e\right )} \,d x } \]

input
integrate(1/sin(f*x+e)/(a+b*sin(f*x+e))^(1/2)/(c+d*sin(f*x+e))^(1/2),x, al 
gorithm="giac")
 
output
integrate(1/(sqrt(b*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)*sin(f*x + e 
)), x)
 
3.1.50.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\csc (e+f x)}{\sqrt {a+b \sin (e+f x)} \sqrt {c+d \sin (e+f x)}} \, dx=\int \frac {1}{\sin \left (e+f\,x\right )\,\sqrt {a+b\,\sin \left (e+f\,x\right )}\,\sqrt {c+d\,\sin \left (e+f\,x\right )}} \,d x \]

input
int(1/(sin(e + f*x)*(a + b*sin(e + f*x))^(1/2)*(c + d*sin(e + f*x))^(1/2)) 
,x)
 
output
int(1/(sin(e + f*x)*(a + b*sin(e + f*x))^(1/2)*(c + d*sin(e + f*x))^(1/2)) 
, x)